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We present a detailed analysis of all archival INTEGRAL data of the accreting X-ray pulsar 
Vela X-1. We extracted lightcurves in several energy bands from 20keV up to 60keV. The 
lightcurves show that the source was found in very active as well as quiet states. During the ac- 
tive states several giant flares were detected. For these states spectra between 5keV and 120keV 
were obtained. The spectra of the active states were found to be significantly softer than those 
from the quiet states. We performed a statistical analysis of the flaring behavior. The resulting 
log-normal distribution of the intensity of Vela X-1 shows that the source spends most of the time 
at an average flux level of 300 mCrab but also that the distribution extends well up to more than 
2.0 Crab. 



7th INTEGRAL Workshop 
September 8-11 2008 
Copenhagen, Denmark 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribntion-NonCommerctal-ShareAlike Licence. 



http: //pos . sissa . it / 



Statistical analysis of Vela X-l 



ttJ o 



; „ 1 nil 


* hi 

' ft 






111* 

- 


_ ' ' ' ' 1 ' ' ' ' 1 


' ' i ' ' ' ' i ' ' ' ' i 

J tv 


V 


' ' ' i ' ' ' ' i 





6 7 

Days in May 2006 



Figure 1: Lightcurve in the 20-60 keV band and hardness-ratio for data from Revs. 433-440, i.e. from 
2006 May 02 to 2006 May 24. For the hardness-ratio the hard band from 40-60 keV and the soft band from 
20-30 keV was used. 



1. Lightcurve and Hardnessratio 

In this work, the accreting HMXB Vela X-l is studied, which has a pulse period of ~283.5 sec 
and an orbital period of ~9.0days. Overall ~ 3.6 Msec of data were analyzed, taken mostly by 
ISGRI aboard INTEGRAL, which is a large amount of data compared to other measurements. Parts 
of the data have been published already (see 4; 6), presenting lightcurves and spectra. A repetition 
of these works was not intended, instead a more statistical ansatz was taken. In that approach 
events were not regarded individually, but all events were binned together in some kind of common 
grid. Lightcurves for the statistical analysis were extracted from ISGRI in three different energy 
bands: 20^10 keV, 40-60 keV, and 20-60 keV. The time-resolution was chosen to be 283.5 sec to 
average each data-point over one pulse period in order to eliminate these fluctuations. Fig. [I] shows 
the lightcurve of the 20-60 keV band for data from the revolutions 433^140 as an example for the 
analyzed data. The vertical dashed lines show the beginning and the end of the eclipse according 
to the ephemeris data from (4), and the dash-dotted line shows the respective center of eclipse. The 
lightcurve was extracted using ii_light. 

The second panel of Fig. [I] shows the hardness-ratio for the respective bins. The hardness-ratio 
is given on the timescale of Science Windows (ScWs) to smooth out the noise. The hardness-ratio 
# is calculated as 

# = — (1.1) 
H + S 

where H is the countrate in the hard band between 40-60 keV and S is the countrate in the soft 
band between 20-40 keV. A larger value of #, i.e., a less negative value, is called harder, where 
as smaller values of # are called softer. The hardness-ratio is a gauge for the physical processes 
taking place and is indicating the shape of the spectra. In Fig. [T] it is evident that in the bright states 
Vela X-l becomes distinctly softer. The two main parts of interest are at May 01 - May 3.5 and 
May 10 - May 12. We note, that in the bright part of the second time range, from May 11 - May 
12, Vela X-l does not show distinct flaring, but is overall almost three times as bright as in first 
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part, from May 10 - May 1 1. To analyze this behavior more closely, we took a look at the spectra 
in these time ranges. As the exposure times are relatively small, we combined the spectra from 
both active states of Fig. [I] and the spectra from both quiet states. The active state is shown in blue, 
while the quiet state is shown in red. The spectra were fitted in XSPEC (1) independently but using 
the same model, namely 

1(E) = C ■ (/phabs ' ^cont) " ^gauabs ' ^gauabs (1-2) 

This model describes a photoabsorbed continuum, superposed by two Gaussian absorption lines, 
one to model the CRSF at ~58.6 keV, the other to model a yet unexplained feature at ~10 keV. The 
constant C is needed, as spectra from two different instruments, namely JEM-X and ISGRI were 
used. The calibration between these instruments is inconsistent, so that slightly different absolute 
fluxes are measured. All other parameters of the model spectrum are the same for both instruments. 
The soft part of the spectra was extracted using JEM-X from 3.0 keV to 26.5 keV, the harder part 
emerges from ISGRI data between 19.0 keV and 130.0 keV. The second and third panel of Fig. || 
show the respective residuals as reduced % 2 value. As shown in Fig. ^| the photo index of the quiet 
state is about 1.5 times as high as the one of the active state. Following (2) this is due to the fact that 
less photons were Comptonized to high energies in the accretion column of the neutron star. This is 
the case, when more material is accreted, e.g. when the absorption coefficient rises, as the photons 
are either emitted or absorbed before they can be scattered to high energies. A larger accretion rate 
additionally leads to an increased brightness in X-rays. This theory can consequently explain the 
effect that the spectrum softens while Vela X- 1 is getting brighter. Another hint for an increased 
absorption coefficient is the fact that in the active state, the Nu value of the photoabsorption part of 
the model is higher, which is a direct effect of more material in the vicinity of the neutron star. The 
energy and the shape of the CRSF does not change within the uncertainties, so that an influence of 
the magnetic field on the shape of the spectrum can be ruled out. 

2. Flaring behavior 

During their observeration of Vela X-l (4) could measure giant flares in the lightcurve of 
Vela X-l, which are much brighter than the usual flares of the object. To understand if such giant 
flares are explainable by the same physical process as the other, common flares in Vela X-l, we 
analyzed the flaring behavior for the whole data statistically. The flares of Vela X-l have typical 
durations of some hours. To get a good coverage of these flares, the chosen temporal resolution 
of the lightcurve of 283.5 sec is well suited. The lightcurve data were binned logarithmically into 
256 countrate bins after eliminating all datapoints measured during the eclipse of Vela X-l. Fig. || 
shows the histogram for the whole 20-60 keV band. The solid line is a Gaussian fitted to the data, 
representing a log-normal distribution. 

Fig. H] shows that the flaring follows closely a log-normal distribution. Even though the system 
is quite different, it is interesting to note that (8) for short time-scales and (5) for long time-scales 
found a similar log-normal distribution when analyzing the flaring behavior of Cygnus X-l. These 
results were confirmed by (3). A log-normal distribution can be achieved by multiplying N random 
subprocesses. The emerging distribution will be log-normal distributed (9). Even if the log-normal 
distribution is a rather good description of the histogram, small offsets at especially high and low 
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Figure 2: Spectra combined of JEM-X and ISGRI data. The blue spectrum shows the active state, the red 
one the quiet state. The second and third panel are the according residuals. 

countrates are apparent in Fig. [|. These offset mean that the process is not produced by completely 
multiplicative components, but small additive components are a not-negligible. To improve the fit, 
two additional Gaussians where added to the model. This improves the fit drastically, as shown in 
Fig. ||. The probability to measure extremely bright flares with more than 300 counts sec -1 in the 
20-60 keV band is with one Gaussian ~0.10%, while when using 3 Gaussians it is ~0. 16%. This 
is a rather small number but far from impossible, so that we can conclude that bright flares are rare 
but not once-in-a-lifetime events. 

A fit with three independent Gaussians was also done for the logarithmic binned lightcurves of 
the soft 20-30 keV band and the hard 40-60 keV band, see figs. |5| and |[ Gaussian #2, describing 
the brighter part of the distribution is even more pronounced in the data of the soft band than it 
is in the data for the overall band. On the other hand, no Gaussian is needed to model the active 
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Figure 4: Histogram of the lightcurve of the 20-60 keV band, binned to 256 bins. The red curve shows the 
best fit consting of 3 Gaussians. The blue curve shows the Gaussian #1, the violet one Gaussian #2 and the 
green one Gaussian #3. 
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Figure 5: Histogram of the lightcurve of the 20-30 keV band, binned to 256 bins. The red curve shows the 
best fit consting of 3 Gaussians. The blue curve shows the Gaussian #1, the violet one Gaussian #2 and the 
green one Gaussian #3. 

end of the histogram for the hard band. Even more, the main Gaussian #1 models the distribution 
almost perfectly, so that the two other become very small. Gaussian #3, responsible for the quiet 
end continuum is almost constant throughout all energy bands, only its overall strength is smaller 
in the hard band. When comparing the peak values of the Gaussians the power-law-like spectrum 
of Vela X- 1 must be regarded, which means that the hard band has intrinsically lower countrates 
than the soft band. Assuming a factor of Soft / Hard of about 5.6 leads to a value of log(5.6) = 1 .73 
which must be added to the center value of the Gaussians of the hard band. This approximation 
would lead to a value of ~ 3.7 for the main Gaussian of the soft band, which is near the actual 
measured value. 

To compare our findings to more measurements data from the All-Sky Monitor (ASM) aboard 
the "Rossi X-Ray Timing Explorer" RXTE was analyzed. The ASM performs 90 min dwells on 
a irregular basis onto Vela X-l, providing a good statistic in the 2-10keV energy band. We per- 
formed the same analysis of the ASM data as of the ISGRI data, with binning the logarithm of the 
countrate to 256 bins after removing datapoints in the eclipse. A model of three individual Gaus- 
sians was then fitted to the data. The 2-10 keV range is clearly softer than the soft band of ISGRI 
so that an even stronger bright excess would be expected. Fig. [7] shows, however, that no bright 
excess is measured. Even more the second Gaussian, responsible for modelling the bright excess 
in the INTEGRAL data has vanished completely. One possible explanation for this behavior would 
be that the photo-absorption in the active state is drastically increased. As the cross-section of the 
photo-absorption is declining fast with increasing energy, X-ray photons of about 20 keV are only 
slightly absorbed, while the 2-10 keV photons are strongly absorbed. The increased absorption is 
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Figure 6: Histogram of the lightcurve of the 40-60 keV band, binned to 256 bins. The red curve shows the 
best fit consting of 3 Gaussians. The blue curve shows the Gaussian #1, the violet one Gaussian #2 and the 
green one Gaussian #3. 



evident in the spectrum shown in Fig. |[ This increased absorption means that even very bright 
flares can not be seen in the soft band, because the soft X-rays can not penetrate the surrounding 
medium. A higher absorption in the active states can easily be explained by the fact that more 
material is accreted onto the neutron star, which means that it moves through denser parts of the 
stellar wind. Another reason that no bright excess is apparent in the histogram of the ASM might 
be the inadequate sampling rate of the ASM. As shown by (7) flares which were clearly visible in 
INTEGRAL data do not show up in the ASM data. Even more, the ASM erroneously detects fluxes 
on the order of the flaring fluxes while Vela X-l was in eclipse (4). This reduces the trust which 
should be put into the histogram of the ASM data further. 

We already noted that brighter parts of the lightcurve are softer. This agrees very nicely with 
the observation that the second Gaussian is vanishing in the bright part of the histogram for the 
hard band, while it is undoubtedly clear in the soft band. This result could mean that two different 
physical processes are at work when producing the X-rays of Vela X-l. One process might be 
responsible for the most of the variability, having the same strength in all energy bands. This 
variability could be due to small variations in the accretion stream, without changing the overall 
configuration of the system. The other process just produces soft X-rays, but with distinctly higher 
countrates. This process is not so common as the first one and might be explained by extremely 
increased accretion which than increases the optical depth in the accretion column dramatically. As 
explained above, this would lead to a softer spectrum. Much more work on this topic is necessary 
to understand the physical processes. 
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Figure 7: Histogram of the lightcurve of the ASM of RXTE in the 2-10keV band, binned to 256 bins. The 
red curve shows the best fit consting of 3 Gaussians. The blue curve shows the Gaussian #1, the violet one 
Gaussian #2 and the green one Gaussian #3. Gaussian # 2 as vanished completely and the bright flank is 
perfectly modelled with only one Gaussian. 
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